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Abstract 

We study the phenomenon of spatiotemporal stochastic resonance (STSR) in a chain of diffusively 
coupled bistable oscillators. In particular, we examine the situation in which the global STSR 
response is controlled by a locally applied signal and reveal a wave front propagation. In order to 
deepen the understanding of the system dynamics, we introduce, on the time scale of STSR, the 
study of the effective statistical renormalization of a generic lattice system. Using this technique we 
provide a new criterion for STSR, and predict and observe numerically a bifurcation-like behaviour 
that reflects the difference between the most probable value of the local quasi-equilibrium density 
and its mean value. Our results, tested with a chain of nonlinear oscillators, appear to possess 
some universal qualities and may stimulate a deeper search for more generic phenomena. 

PACS numbers: 02.50.-r,05.40.-a,05.50.+q,87.10.+c 
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Since the appearance of [l] the phenomenon of stochastic resonance (SR) has become a 
popular field of research. A great deal of experimental as well as theoretical and mathemati- 
cal work has been devoted to the study of the phenomenon in different systems (for reviews, 
see j^l I2I 0, There has been a particular emphasis of its relevance and importance in 
biology and medicine [f| where noise in general, and SR in particular, play a surprisingly 
constructive role 131 . The abilitv of SR to generate order from disorder IS 



especially of relevance in the context of pattern formation mechanisms that are enhanced 
by noise ^|. The discove ry o f an enhancement of the effect by the coupling of nonlinear 



>veryo: 

m 



oscillators into an array |H| has brought new insight to studies of SR. This effect in a wide 
sense is known as spatiotemporal SR (STSR) 0|. An explanation of this effect has previously 
been described as "collective spatiotemporal motion" and "optimal spatiotemporal synchro- 
nization". In spite of much progress, the precise kinetic details of such a synchronization 
remain without an appropriate study. 

We conjecture that among the kinetic details of STSR, wave front propagation plays a 
prominent role ( H, Q] are good background references in this regard). We investigate 
the related problem of inducing and controlling global spatiotemporal order in a chain of 
diffusively-coupled bistable oscillators by a locally applied signal. The ability to induce STSR 
throughout the chain by applying a local signal to a small part of the chain is precisely shown. 
This can be regarded as an element of a spatial signal transmission, and possibly gives a 
new design freedom to modelling biological and biomedical problems. We also investi gat e 
the effective statistical renormalization of the steady states of a generic lattice system [24]. 
This renormalization reflects the difference between the most probable and mean values of 
the local quasi-equilibrium density which is a result of time averaging [25]. This leads to 
a new observation that the system, on the time scale of STSR, exhibits a bifurcation-like 
behaviour. It also gives a criterion for the noise intensity depending on the coupling in the 
chain. Both our results appear to have a certain universal quality and may stimulate a 
deeper search for generic phenomena, e.g. in a chain of FitzHugh-Nagumo equations jlfil ]. 

Consider a chain of overdamped oscillators with diffusive coupling of constant K > 
and a bistable on-site potential V (y) = —my 2 / 2 + y 4 /4. We assume that the system is 
influenced both by external random noise of intensity D, which involves a set of independent 
generalized Gaussian random processes {£ n (^)} with two characteristic cumulants, — 
and (£m(t)£n(t')) = 5mn<5 (t — t'), and a deterministic signal S(t) applied locally to a part 
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of the chain, Sn (t) = {s(t) if n G and otherwise}. In what follows we fix 

the particular and simplest form of the external signal, s(t) = Acos(W), together with 
dimensionless amplitude A = 0.025, and frequency uj = 5tt • 10~ 5 that actually set one of 
the timescales, T s = 2ttuj~ 1 = 4 • 10 4 . The other characteristic timescales are: the relaxation 
time of the chain to Gaussian fluctuations in the vicinity of one of its stable steady states, 
T r (in our case T r ~ (2m) -1 ~ 1); the waiting time of the initial birth of an "instanton", 
Tk (Kramers' time) j9( - here we do not explicitly consider this timescale but, fortunately, 
there are in-depth studies of the problem Q, (a posteriori one can say that Tk is 
considerably shorter than the principal timescale here) and also the timescale related to any 
wave-front propagation in the chain, T w . 

The corresponding chain stochastic differential equation (SDE) in dimensionless variables 
has the form, 

y n = KAy n - V (y n ) + V2Di n {t) + (t) , (1) 

with n on the chain N in the integer lattice Z, n G N = {1, . . . , k+1, . . . , fc+M, . . . , N}, = 
{k + 1, k + 2, . . . , k + M} , ^ M ^ N; and Ayi = y i+ i — 2y i + y i _ 1 . Two different topologies 
are possible for the chain: either the ends are connected or not connected. We take the latter 
case with Neumann boundary condition. 

The corresponding physical context of (0) is in the Smoluchowski kinetics of a harmon- 
ically coupled chain of particles with transverse displacements {y n }- The underlying free 
energy functional and the corresponding deterministic part of the dynamics, without sig- 
nal, are: T (y) = J2( n ) \ K iVn - Vn-\f + V (y n ) , y n = -d? '/dy n . This interpretation can 
considerably facilitate the understanding of the dynamic behaviour of ([Tjl. Especially since 
even diffusively-coupled, chain oscillators a ppe ar surprisingly difficult to analyse consistently 
from a rigorous mathematical perspective [18j | . 

We study a spatially discrete model because chains or array structures are often of rel- 
evance in biology. Besides biology there are also crystal lattices |l7(. Moreover interesting 
dynamical effects exist in discrete models that are not present in their continuous analogs - 
e.g. the propagation failure of travelling waves [3], and breathers [3]. Even in population 
dynamics it was recently demonstrated that important effects due to the discrete nature of 
organisms may be entirely missed by continuous models 2o| . 

To introduce induced STSR phenomena, we first present a representative numerical sim- 
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ulation of (0) that is performed longer than all the characteristic timescales. We use time 
steps 0.01 and 0.001 that are considerably shorter then the principal timescale. Further 
shortening of the step does not change the result. As shown in Fig. [T] the system can indeed 
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FIG. 1: Locally induced global spatiotemporal pattern in the chain of bistable oscillators; param- 
eters are specified in Fig. |^1 



represent the well-recognized phenomenon with a local signal applied to only 1/6 part of the 
chain. The collective variables, Y = N -1 Yln=i Vn an d Y' = M _1 Yl^k+i Vn, also adequately 
and legibly reflect the key features of the effect as shown in Fig. |21 
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FIG. 2: (a) Dynamics of the collective variables Y (solid) and Y' (dash) with parameters N = 300, 
M = 50, k = 125, K = 15, D = 0.1, and m = 0.25; (b) Fine structure of the transition region in 
(a) reveals a wave front propagation. 



The generic features of STSR, i.e. those related to a diffusively coupled chain of bistable 
oscillators, have reasonable prototypes going back as far as the pioneering paper [jj. The 
main result from [l| is that, for transition times, only the energy of the unstable instanton- 
like spatially inhomogeneous steady state solutions, and not the total energy barrier for the 
chain, is of importance cf. [21] . However, in fact, there has been no discussion at all of the 
transition kinetics, which appear after the birth of the "instantons" . 

The transition kinetics are most likely related to, but different from, another interesting 

n 

phenomenon - travelling waves 18]. It is evident from the physical interpretation of (JTJ, that 
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a travelling wave solution can appear only if the symmetry of the bistable potential is broken 
and there is an energy gap between the two stable states, and the energy flux is compensated 
for by dissipative forces. This is not the case for the deterministic part of (0) with A = 0. 
However, if A ^ then the external signal periodically breaks the symmetry and creates 
an absolute minimum at one of the wells of the underlying potential, and a travelling wave 
front can develop on a sufficiently long timescale T s . The situation changes favourably in 
the case of a chain because the translational symmetry is broken. With the chain, starting 
from the physical interpretation above, we can imagine an unstable instanton-like steady 
state solution p. Its characteristic lifetime introduces a new time scale that is required to 
be shorter than T s for the continued existence of the induced STSR effect. Our numerical 
simulations justify a posteriori that this condition is perfectly realizable. 

Once the underlying kinetic mechanisms of STSR are understood, we can explore other 
interesting features. Here we focus on one that has been overlooked in previous studies, viz. 
the shift of the stable spatially homogeneous steady states. 

Let us consider the specific lattice SDE 

Vi = KA Vl - V (y<) + V2D& (t) , i G Z, (2) 

and evolve the optimal Gaussian representation of this equation (extending to the case 
of a lattice system as a starting point for further analysis; alternatively one can start with 
an averaging principle j^j]), 

ili = KA Vi - [ ai (t) + k (t) 5yi) + V2D& (t) , i G Z, (3) 

where and 6j are obtained by a minimization procedure in respect of the mean-square 
error functional J = ([V (yi) — (a« + biSyt)} 2 ') for all i G Z, where 8yi = yi — (yi). The 
representation Q of © is especially appropriate since the ex ante fluctuations are almost 
Gaussian near the stable steady states and typical exit paths of (0). On the timescale 
T s the internal deterministic and random oscillations are very fast and can be considered 
adiabatically following with the external signal. 

To actually obtain and 6j, which is an intractable problem since it requires the exact 
solution of 0, we consider a self-consistent approximation scheme combining the minimiza- 
tion of J together with the solution of Q. Thus we replace averaging according to © by 
its Gaussian approximation according to (J3J), y — (y) Gaussian = (y)- As a result, we obtain 
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the self-consistent set 
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(4) 
(5) 



Further, using the Furutsu-Novikov formula [231, and references] we explicitly obtain Oj = 
(3£ 2 (z/i) -m)V% + Vh h i = (3&2{yi) - m) + 3f/f, where & 2 {yi) = (yf) ~ vl is the second 
cumulant. Thus (HJ) takes the form §i = KAiji + (m — 3^2(2/1)) Vi — yf, involving the effective 
potential function V eS , V^ s (&) = - (m - 3jft 2 (yi)) ■ Since # 2 (j/i) > 0, V cS > V, always. 

Lastly we consider the principal problem of an explicit calculation of ^{Vi) and solve it 
under certain hypotheses. Consider the spatial correlations, K mn (t) = (Sy m (t)8y n (t)) and 
suppose that the fluctuations tend to their steady state via a stage of spatial homogenization 
with the ansatz K mn (t) = K m - n (t) = K r (t). The known equilibrium solution of (0) with 
probability density Poa(y) exp[— J 7 (y)/D] has the property of spatial homogeneity, but: 
Does the property still persists as time evolves? To facilitate understanding, consider the 
process of the formation of 

Poo ^ 

result of time averaging, and fix two limiting cases 
related to low and high levels of the noise intensity. As D I the system spends most 
time in a potential trough, occasionally passing from one to the other. Being in a trough it 
has time to form a local quasi-equilibrium density that reflects the local asymmetry of the 
underline potential. As t — > +00, a sum of the local densities is formed to satisfy the global 
symmetry condition. Under a high level of noise, the rate of passage from one trough to the 
other is frequent enough in order to rapidly form the mean value y = 0. The qualitative 
picture described above is linked to the particular time scale of the problem in question. 
Favourably for the ansatz, we are interested in averaging over a time scale about T s that 
becomes apparent in the local quasi-equilibrium density. Test numerical simulations are also 
essential in order to reinforce the ansatz. 

Using (JSJ together with the Furutsu-Novikov formula, we obtain a dynamical equation 
for K r (£): k r = 2KAn r — 2b(t)n r + 2D5 r Q, r G Z, - and as a result, the equation 
for the steady state correlation function, KAn r — bn r + D5 r o = 0, r e Z , with the 
natural asymptotic conditions lim^-too K r — 0; b is still unknown. Substituting n r = 
> 1, we obtain the set of algebraic equations, corresponding to r = and r > 1 : 



[2K(1 - t- 1 ) + b]A = D, t 2 - 2(1 + b/2K)t +1 = 0. The last equation has two different 
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FIG. 3: Spatially-homogeneous steady states y ^ as points of intersection of a characteristic 
surface y = y(D,K) and two coordinate planes of the control parameters {D = 0.1 and K = 15 
are marked). Instead of the two original steady states ±0.5 there are now four for specific values 
of D and K, and zero otherwise. 



roots, t± = (2K)' 1 (2K + b± \/b(AK + b)) , connected by the relation, t+ ■ t. 
t = t+ > 1, this means that t- = t^ 1 < 1. Therefore finally we obtain 

|r| 

D 2K +b + \/b{4K + b) 

K r 



1. Since 



y/b(AK + b) 



2K 



In particular, Kq — D/ ^Jb(AK + b). Further, we can use k in combination with (y) = 
to characterize the effective spatially-homogeneous steady states. Excluding b (= 3k,q — m + 
3y 2 ) from this set, we arrive at two cases: 



(a) y = 

(b) f 



0. 



(3kq — m) (AK + 3kq — m) = 



2 • 



m — 3kq, 



D 2 

(m — 3k ) (2K + m — 3k ) = 

4:Kn 



Observe first that while (a) always gives y = and k ^ m/3, the set (b) has no real 
roots y («o > 0) for a certain range of values of D and K (see Fig. |2J). In other words, 
the stable steady states disappear and only y = remains in this range. We could claim 
that this occurs with a bifurcation-like behaviour. To identify this observation, we carry out 
numerical experiments on the system (Q) varying the noise intensity D over a wide range. 
The response to the external signal provides evidence of the effect as shown in Fig. 01 which 
is enhanced further with larger N. In the case (c) the stochastic resonance pattern, adjusted 
for the renormalization, is recognizable, but in the case (a) only simple oscillations with the 
frequency of the external signal are visible around y = 0. 

To further understanding of this observation, one can consider the equilibrium probability 
density of (0). As D j 0, Poo(y) is concentrated at two spatially- homogeneous absolute 
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FIG. 4: Three output signals (collective variable Y, N = M = 500 ) with (a) D = 0.5; (b) D = 0.24; 
and (c) D = 0.1; correspondingly, above, near, and below the characteristic surface y = y(D,K) 
at K = 15 (see Fig. |3J); effective stable steady states: (a) 0, and (c) ±0.43. The difference appears 
to be not only quantitative but also, more importantly, qualitative. 

minima of the potential and forms the most improbable configurations around unstable 
steady state y = 0. The homogeneous equilibrium mean value (y)oo = is fixed by the 
invariance of Poo{y) under the transformation y — > —y, and eventually coincides with the 
unstable steady state. By increasing D we can change the sharpness of the equilibrium 
density profile, but not the topology. It should be observed here that Poo(y) represents the 
average over all sample paths of (J2J). If we consider a single sample path then the local density 
profile of (J2J) is formed by time averaging, initially in the vicinity of an absolute minimum 
of the potential reflecting the local potential asymmetry. After a sufficiently long period 
the system with probability 1 will pass through the low-probability range to the vicinity of 
the other absolute minimum of the potential, and so on. The rate of this process rapidly 



increases with D 2l|. The rate of development of a difference between the most probable 
and mean values of the local density is also increasing with D since, within the same time 
interval, more different states are realizable. It is remarkable that the time interval related 
to the signal as well as the inherent nature of SR is favourable to observe the mean values 
or, in other words, to feel the effective potential. The collective variable Y, which is a sort 
of site average, allows us to visualize this effect if the chain is long enough. 

Finally, if we set D = 0.1 and \fm = 0.5, we obtain y m ±0.43. This is close to the 
observed value obtained by direct simulation of (0) (with N = M = 300 and 500). The 
agreement improves with larger values of N. There also exists a pair of steady states close 
to y = 0: y « ±0.1, - but these are not clearly identified in the simulation of STSR. 
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